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We study the collective modes of the quantum Hall states in undoped bilayer graphene in a 
strong perpendicular magnetic and electric field. Both for the well-known ferromagnetic state that 
is relevant for small electric field E± and the analogous layer polarized one suitable for large E± , the 
low-energy physics is dominated by magnetoexcitons with zero angular momentum that are even 
combinations of excitons that conserve Landau orbit als. We identify a long wave length instability 
in both states, and argue that there is an intermediate range of the electric field < E± < 

(2) 

Ec where a gap less phase interpolates between the incompressible quantum Hall states. The 
experimental relevance of this crossover via a gap less state is discussed. 

PACS numbers: 71.35.Ji, 71.70.Di, 73.43.Lp, 75.30.Ds 



I. INTRODUCTION 



Recent magnetotransport experiments^ - — on bilayer 
graphene 5 crystal in Bernal stacking (Fig. [1]) have shown 
that its characteristic, eight-fold quasi-degenerate (spin 
a =t, |, valley £ = ±1, and Landau orbital n = 0, 1) 
central Landau band (CLB) is split at all integer values 
of the filling factor v = ph/eB±. Quantum Hall ferro- 
magnetic states have been suggested 6 to explain these 
states. At v = 0, on the other hand, the longitudi- 
nal resistivity p xx diverges beyond a sample-dependent 
threshold value of the magnetic field both for decreasing 
temperature and increasing field. Similar behavior has 
been observed in monolayers^ and in particle-hole sym- 
metric semiconductor systems. 8 While this is unusual, it 
does not rule out quantum Hall physics, because Laugh- 
lin's gauge argument 9 connects a vanishing longitudinal 
conductivity a xx to the quantized Hall conductivity a xy , 
which is consistent 10 with a divergent p xx = p yy for 



P: 



xy 



0. While Zhao et 



- found that the gap at v = 
B±, as expected of Coulomb 



depends on the field as 

interaction effects, Feldman et aZ.,i using high- mobility 
suspended samples, measured a gap that opens linearly 
with B± and hardly depends on Bn. The latter finding 
suggests that many-body effects dominate the Zeeman 
splitting. The observed linear B± dependence of the gap 
up to a rather high threshold value has been explained 
by Nandkishore and Levitov-ii and Gorbar et ali^ taking 
dynamical screening into account. Thermally activated 
transport across the ex B± gap then explain o 11 ^ 2 the ex- 
ponential growth of resistivity p xx in B^/ksT found in 
Ref. H. 

At zero magnetic field, recent experiments^ have con- 
firmed the emergence of a band gap 14 if a perpendicular 
electric field is applied. In the presence of a perpendicu- 
lar magnetic field B±, the band gap affects the Landau 
levels (LL's); but as exchange energy considerations are 
fundamental in the quantum Hall regime, the low-energy 
physics at integer filling factors is determined by the in- 




FIG. 1. (Color online) The crystalline structure of bilayer 
graphene in Bernal stacking. A, B (A, B) denote inequivalent 
sites on two hexagonal carbon sublattices on the bottom (top) 
layer. Hopping amplitudes are indicated. 



terplay of the electron-electron interaction, the Zeeman 
energy Az = gpBB±, and the inter layer potential en- 
ergy difference Al = edE± — An/(2eo) between the lay- 
ers (An is the electron density imbalance, d = 0.335 nm 
is the distance between the layers, g is the gyromagnetic 
factor, and ps is the Bohr magneton). In a mean- field 
approximation, the ground-state of this system sponta- 
neously break o 12 i 14 the (approximate) spin and valley 
symmetry^ Obeying variants of Hund's rule, the ground- 
state is ferromagnetic in the A^ <C Az limit and valley- 
polarized in the Az, ^> Az limit. As the Coulomb 
interaction, however, is quite strong in comparison to 
the LL splitting, [e 2 / / (47re e r ^)] / ' E 2 w 11.7/[^ (T)e r ] 
(E2 is the greatest Landau level energy difference, c.f. 
Sec. HI1 below), screening by inter-LL transitions might 
be important . n i 12 (Here £3 = yh/eB± is the mag- 
netic length, and e r is the relative dielectric constant of 
the bilayer graphene sheet and its immediate environ- 
ment.) Thus, when studying the excitations of bilayer 
graphene on the a xy = plateau, we may start from 
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these symmetry-breaking quantum Hall states. We shall 
regard, however, the perpendicular electric field E± in- 
stead of Al as the tunable parameter. 

The paper is organized as follows. In Sec. [Ill we review 
the Landau- level structure of bilayers and recall the sym- 
metry breaking ground-states that are relevant at charge 
neutrality. In Sec. IIII1 we define the type of excitations 
we study. In Sec. IIV[ our results are presented, and in 
Sec. El their significance is discussed. Technical details 
are included in the Appendixes. 



II. LANDAU LEVELS AND GROUND-STATES 

In the vicinity of the valley centers corresponding to 
the K (£ = 1) and K' (£ = —1) first Brillouin zone 
corners, the electronic structure of bilayer graphene is 
well described by the tight-binding Hamiltonian 15 
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where 7r = p x -\-ip y and p = -?W-eA, v = \ /r 3ajo/2h « 
10 6 m/s is the intra-layer velocity, and 71 « 0.39 eV is the 
inter-layer hopping amplitude. This Hamiltonian acts in 
the basis of sublattice Bloch states [ipA, ip§, ^] m 
valley K and [ip^^A^B^x\ m va ^ e y K' • Having a 
suspended sample in mind, we will assume e r = l,g = 2. 

We emphasize that our main results follows from the 
four-band model of Eq. (pQ); the low-energy two-band 
model 14 is only occassionally referred to for contrast. 
We have neglected 73, 74, A 7 of the Slonczewski-Weiss- 
McClure model 16 as usual. In particular, Ref. [l5| has 
shown that 73 ~ O.I70 has no significant effect on the 
Landau levels of bilayer graphene, which suggests that 
orbital mixing due to 73 is negligible. Using typical val- 
ues from the literature, 74 « 40 meV and A' « 7 meV 
are small in comparison to the Coulomb energy scale 
e 2 /(47re e r £ B ) = 56 meVy/B[T]/e r . 

Using the gauge A = —BySt, the Landau levels and 
orbitals are obtained from the ansat z 17 i 18 



*p0q = [^0g,0,0,0] T , 



(2) 



iftnaq — [Aa^Vnq, ^a^Vn-2,q, iC^TJ n -i iq , iD^TJ n -i iq \ , 

where n > 2; 1 < a < 4 (E na < E n ^ a +i), and r] nq are 
the single-particle states in the ordinary two-dimensional 
electron gas with quadratic carrier dispersion, 

iqx-(y/£ B -q£B) 2 /2 / \ 

Vnq (r) = - ^H n ( -j- q£ B J , (3) 

^2n^2 n n\£ B Vb J 

and H n is a Hermite-polynomial. Introducing the dimen- 
sionless quantities 

t = ji£ B /vh, S = A L £ B /2vh, 



the n — 1 Landau levels Ei a £ = hvei a ^/£ B ^ a = 1,2,3, 
are obtained from the ensuing secular equation, and the 
orbitals are determined by 
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where N\ a is the appropriate normalization factor. The 
n > 1 Landau levels E na ^ = hve na ^/ '£b, a = 1,2,3,4 
follow similarly, The orbitals are specified by 
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D 



^2n - 2 



2n 



2nN n 



-a 



with the appropriate normalization factor N na . 

For the central Landau bands, in particular, we obtain 



E- 



l,a=2,£<7 
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States tpo^aq and ipi, a =2,£aq form the central Landau 
band (CLB) octet. For n = (n = 1) the states are 
completely (predominantly) located on A sites in valley 
K and B sites in valley K' , making valley equivalent to 
layer in the central Landau band. 

With this notation, the ferromagnetic state is£ |f.m.) = 

^o,£,tg^!,cK=2,£,tJO)' anc ^ ^ e valley (layer) polarized 

one is |v.p.) = n ? a a o^=i,^ a i,a=2^=i,aal°)- Here 1°) 
is the "vacuum" of an infinite number of filled valence 
Landau bands (n > 2, a = 2, 3, all combinations of £ = 

±1 and a =t,40, ano ^ ^ia^crp = ^Ivp creates a particle in 
state ipNp- Al differs in the two states: 



A f L m - = eE ± d, 



A L ^ = ed(E ± -E ), (4) 



where E = e 2 d/ (7re e r £ 2 B ) « 8.75 x 10 6 £[T]/e r V/m. 
Thus in state |v.p.) A^ is reduced from its unscreened 
value by ~ 34B[T]/e r K, which is greater than the Zee- 
man energy Az ~ 1.4S[T] K. 



III. MAGNETOEXCITONS 

The low-energy excitations of quantum Hall states are 
magnetoexcitons, which are obtained by promoting an 
electron from a filled Landau band to an empty band 19 . 
These neutral excitations have a well-defined center-of- 
mass momentum Q. They are collective, but also in- 
clude widely separated particle-hole pairs in the Q — >• oo 
limit. The latter limit determines the transport gap un- 
less skyrmions formSS. Magnetoexcitons are created from 



the ground-state by operators^ ^jvjv'(Q) 



*U'(Q) 



In the low magnetoexciton density limit, the interac- 
tion between magnetoexcitons is neglected; then it is suf- 
ficient to diagonalize the mean-field Hamiltonian matrix 



<$(Q) 



^ m ,(Q)V^ NN ,(Q)\0) 



(0\V\0), (6) 



where N = (n, a, £, <j) [N f = (n', a', a')] specifies the 
Landau band where the particle (hole) is created; 5nn' = 
Saa'Szz'Snn'tiaiai'- Taking the four-spinor structure of the 
Landau orbitals is the only nonstandard step. Details are 
delegated to the Appendixes. 

Magnetoexcitons carry spin and pseudospin (valley) 
quantum numbers, as defined by the particle and hole 
Landau bands involved. In addition, an "angular mo- 
mentum quantum number" l z = \n\ — \n'\ can be defined 
as usua h 19 i 22 This quantity, however, is exactly conserved 
only in the Q —> limit. 




IV. RESULTS 
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FIG. 2. (Color online) Excitations of the neutral graphene 
bilayer in perpendicular fields in the absence of an inter- 
layer energy difference A l • Even and odd refer to the modes 

i^'=o' eVen (k) and i£/^Q' odd (k), which at low momenta tend 

, p p 

to the even and odd linear combinations, oc ^ (o^uo^'t) ± 
,p p 

^ (l'^i^t)' °^ exc itations that conserve Landau orbitals, re- 
spectively. The inset illustrates qualitatively the band split- 
tings by Al for unrealistic Zeeman energy (B « 10 4 T) for 
visibility. 



The magnetoexcitons of the ferromagnetic state |f.m.) 
are characterized by their pseudospin P, P z , and angular 



FIG. 3. (Color online) Excitations of the neutral graphene 
bilayer in perpendicular fields at various parameters B± 
and Al/Az as indicated on the panels. The energy unit 
e 2 /(47reoe r ^s) already scales with B; the additional change 
shown here is due to the B dependence of the Landau or- 
bitals. 



momentum l z quantum numbers. The latter is l z = 1 

when an electron is promoted from an n = orbital to 

n = 1, and l z = — 1 in the converse transition. (For |v.p.), 

spin and pseudospin are interchanged.) The spectra are 

shown in Fig. [2] for Al = 0. The l z = excitations have 

two branches, E^Q even (k) and E^Q° dd (k), each with 

a pseudospin triplet and a pseudospin singlet band. In 

the long- wavelength limit, the i^ p ^Q' even (k) mode tends 

to ~ Az + P z Al (apart from a small P^-dependent cor- 

+p p 

rection). This excitation is ex \I> (o&uo£'t) 



w (ia,i£'t) 

in this limit, while the gapped Ep^Q° dd (k) mode is 



±1 branches are 



^ *VaV't) * (i'a,iet)- The ^ 

gapped and degenerate, in compliance with the global 
particle-hole symmetry of the system. 

For nonzero A^, each band is split according to its 
pseudospin and spin quantum numbers as sketched in 
the inset^ of Fig. [2l For generic Al, the dispersion of 
magnetoexcitons also depends directly on Bj_ (not just 
through £b) because of the 5 = AlIb/^vH dependence 
of the orbitals; c.f. Fig. [3l The excitation energies are 
reduced by roughly 1% to 16% from the two-band value 
| Y / ^e 2 /(47reoe r £ J e), roughly linearly, in the range 1 T 
< B± < 40 T. The density-of-states peak due to the 
flat region in the upper l z = branch and the l z = ±1 
branches might be observable in Raman scattering^ The 
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lower l z = branch is less curved for increasing B±; 
but as Fig. 2] demonstrates, this reduced spin stiffness 
is insufficient to make the charge gap skyrmionic unless, 
perhaps, for huge fields B± > 30 
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FIG. 4. (Color online) Charge gap from widely separated par- 
ticle hole pairs (excitons) and skyrmion-antiskyrmion pairs. 
The latter is independent of the potential difference between 
layers A l ; the former is shown for two different values of Al . 

Starting in the A^ <C Az limit in the ferromagnetic 
state, there is a long-wavelength instability due to the 
pseudospin triplet mode, 

(^(0£=-i,t,oe=U) + ^(i£=-i,t,i£'=H)) ( 7 ) 

at Az ~ A^, where the excitation energy of the mode 
(see Fig. [2]) reaches zero. This is a signal of a transition to 
the other state, which becomes energetically favored at 
this point [the strict equality in Eq. (O holds at k = 0]. 
Figure [5] shows the regions of stability in terms of B± and 
Al/ Az- Conversely, starting from the valley-polarized 
state in the (A^ ^> Az) limit, there is an instability due 
to the mode 

<i A= - 1,evm (k) - 

(^(o£=i,+,oe=-i,t) + *(u=i,^i^=-it)) ( 8 ) 
For realist ics fields, the phase boundary is given by 

(Fig. O 26 

^ « 1 + 0.008755 [T]. (9) 
Az 



V. DISCUSSION 

In experiments, the energy difference between the lay- 
ers, Al, is determined indirectly via the electric field 
E ± . By Eq. © there is a region E^ < E± < E^ 2 \ 
with e[ 1] = Az/ed « 3.5 x 10 5 £[T] V/m and e[ 2) = 




B(T) 



FIG. 5. (Color online) Zero temperature phase diagram of 
undoped bilayer graphene. The red shaded region above the 
thick line corresponds to a compressible state. The contours 
in the ferromagnetic phase show the long- wavelength limit of 
the magnetoexciton gap. 

Az/ed + e/(7T€ £ 2 B ) « 9.1 x 10 6 B[T] V/m, where the 
ferromagnetic state is destabilized by the lowest l z = 
mode, while the valley-polarized state is not yet stable (it 
would overscreen the external electric field and Al > Az 

(2) 

would not hold). For E± > Ec , the valley-polarized 
state becomes the ground-state. 

As |v.p.) and |f.m.) have extremal total Zeeman and 
layer polarization energy, there is no energetically more 
favorable quantum Hall ferromagnet in the i^ 1 ^ < E± < 
E^ range. When E± reaches ^c 1 ^, the ground-state en- 
ergy can be lowered by creating magnetoexcitons. Two 
scenarios are consistent with this observation: (i) a first- 
order phase transition between the ferromagnetic and the 
valley-polarized quantum Hall states, and (ii) a crossover 
via compressible states in the E^ < E± < E^ range. 
At zero temperature, this issue is decided by the energy 
comparison of the Maxwell construction that interpolates 
between the energies of the two states |f.m.) and |v.p.), 
and the state that is obtained from these quantum Hall 
states by populating the k = state of magnetoexcitons 
of the critical mode of Eq. 0. Formally, the latter can be 
conceived as a Bose-Einstein condensation of magnetoex- 
citons, but it merely amounts to a gradual change of the 
many-body ground-state in the crossover, with no quasi- 
particles keeping their identity while condensed. If the 
magnon-magnon interaction is neglected, these energies 
are degenerate. If, on the other hand, the interaction be- 
tween magnetoexcitons assumes a van der Waals profile, 
as is typical for the inter-exciton interaction in analo- 
gous systems at vanishing magnetic field ^21 the energy of 
the compressible crossover state is lower than that of the 
phase-separated quantum Hall systems, thus, at least in 
the neighborhood of i^ 1 ^ and E^ where the two-magnon 
interaction dominates over many-magnon effective inter- 
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action terms, the gapless state prevails. The calculation 
of the magnon-magnon interaction is delegated to future 
work. 

Most experiments have been performed^ with a single 
backgate. The gate voltage to counter possible extrinsic 
doping was \V D \ < 0.5 V in Ref. 1 and V D « 4 V in Ref 
The electric field then must be a few volts on about 300 
nm, i.e., - 3 x 10 6 V/m. With B ~ 10 T, this must be 
within the ferromagnetic region. Moreover, the ^-factor 
might be significantly enhanced^. Thus we expect that 
only dual-gated bilayer graphene^^ can be tuned into 
the gapless state. 

More recently, dual-gated two-terminal magnetotrans- 
port measurements were reported by Wietz et al£. A 
transition exhibiting increased conductance was observed 
between two a w plateaus in the two-terminal conduc- 
tance as a function of the external electric field for not too 
small magnetic fields. The threshold field E± of the tran- 

(2) 

sition is about 40% greater than our Ec , and its slope 
as a function of B±, 11 mV/(nm-T), compares well with 
our prediction of 9 mV/(nm-T). Admittedly, the calcu- 
lation we present would imply a wider conducting region 
than that found in Ref. y. Four-terminal experiments 
by Kim, Lee, and Tutuc4 also found a linear dependence 
of the critical electric field with a slightly higher slope, 
12 to 18 mV/(nm-T), and a more stable ferromagnetic 
state. Inter-Landau level screening, neglected in this pa- 
per, may stabilize the ferromagnetic state, in accordance 
with Refs. 11 and HI A more conclusive comparison be- 
tween theory and experiment requires further progress. 

To summarize, we have presented the excitation spec- 
tra of undoped bilayer graphene in the quantum Hall 
regime when the interlayer bias is introduced by a per- 
pendicular electric field. The collective modes of zero an- 
gular momentum that are an even combination of transi- 
tions that conserve the Landau orbital give rise to a long- 
wavelength instability, which precipitates a compressible 
region between two quantum Hall states for a range of 
the perpendicular electric field. 
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Eq. (j6j), where N = (n,a,£,a) specifies a Landau band. 
Suppressing spin and valley labels for simplicity, the field 
operators for bilayer graphene Landau orbitals are 

*(r)= [A n ai Vn q (r),B^Vn-2, q (r), 

nct^aq 

iC a£Vn-l,q(r), iD^T] n - ljq (r)] T CLnataq, (Al) 

cf. Eqs. § and 0, and A^B^C^D^ are set to 
zero whenever redundant for n = 0,1. (That is, we set 
A° ai = for a = 2,3,4, = = D° ? = B^ = 
for a = 1,2,3,4, and A\ = B\ = C\ = D\ = 0.) The 
Coulomb Hamiltonian is 

H=^j dxdx'* t (x)(g)* t (x , )V'(x-x')*(x , )®*(x), 

(A2) 

where V has a matrix structure as it is sandwiched be- 
tween spinors. If f3 runs over spinorial indices (/? = 
1,2,3,4 for sites A, B, A, B in valley K and B,A,B,A 
in valley K' ', respectively), V has the form 



e 2 5 ^5 ^5^5^ 
47reoe r 



£l(£i,/3i),L(6,/32) , h{^ 1 ,(3 1 ),l-L^ 2 ,f3 2 ) 



x - X' 



Vlx-x^ + d 2 



(A3) 



where d = 0.335 nm is the distance between the layers, 
and the layer index L(£, f$) is defined as 

0, if f = 1 and /3 = 1,4, 
P)= { or £ = -1 and /3 = 2, 3, 

1, otherwise. 

Substituting Eqs. (fAlj) and (fA3]) into Eq. (fA2]) yields the 
following terms: 

(i) The single-particle energy difference 

^ N N' N' i^na^a E n > ^> ) , 



This work was supported by the Lancaster University- 
EPSRC Portfolio Partnership. C. T. was supported by 
Science, Please! Innovative Research Teams, SROP- 
4.2.2/08/1/2008-0011. We thank Judit Sari for carefully 
reading the manuscript. 



Appendix A: Mean-field theory of magnetoexcitons 
in bilayer graphene 

While the mean-field theory of magnetoexcitons is 
standard material, we present its adaptation to bilayer 
graphene. The mean-field approach amounts to the di- 

agonalization of the Hamiltonian matrix iJ^^(Q) in 



where 5 NN > = <W%'<W<W . 
(ii)The exchange energy difference 



a ex 



M filled 



««,WC(P)) (A4) 



of the promoted particle in the two bands with all par- 
ticles of like spin and valley. (Ij^InI (p) w ^ be defined 
below.) 

(iii) The dynamical interaction 



w (NN f ) (n \ 

(NN') v^w 



dq 

W) 2 



3 iz-(qxQ) jN N' / \ 
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(iv) A self-energy term 



if ri > n, else F nn /(q) = F*, n (-q). is an as- 

sociated Laguerre polynomial. Notice that E^^ f ^(Q) oc 
^n'-n.n'-n is n °t necessarily diagonal as a matrix indexed 
by Landau level pairs. 



that is associated with the recombination and the recre- 
ation of the particle-hole pair, as traditionally obtained 
from the random phase approximation 19 . Recombina- 
tion, however, is precluded in the studied problem be- 
cause all excitations flip spin in state |f.m.) and valley in 
|v.p.>- 

Here I^nI (p) * s ^ ne sum °f a same layer contribution 
SS(p) and an interlayer one r^(p): 



oN ± N 2( n 

nzNa vp; 



27T6 2 

ep 



azazazazf Ni nAp)f* N2N3 (p)- 



AZC%C2lAZF NlN4 (p)F^ 2 _ 1>N3 _ 1 (p)+ 
C n a \ Al\ AllCZ F^-i^-i (P)F* N2N3 ( P )+ 

Ca 1 1 CS 2 2 C2 3 3 CS 4 4 ^iv 1 -i,iv 4 -i(p)^ 2 _ 1)N 3-i(p)+ 
B2 1 1 ^ 2 2 B2|^:^ 1 _ 2 ^ 4 _ 2 (p)F^ 2 _ 2)N3 _ 2 (p)+ 

B2 1 1 ^^^F Nl _ 2 , J v 4 _ 2 (p)F^_ 1)W3 _ 1 (p)+ 
DZBZBZDZF Nl -i, Ni -i{p)F* N2 - 2 , N3 - 2 {p)+ 
DZDZDZDZF^^^F^^ip)], 

2ire 2 



T NiN 2 1 \ 
J-NzN* \P) 



-pd 



ep 



[Al\BllBZAl\F NlN MF* N2 - 2 , N3 - 2 (p)+ 

Bl\ A Z A Z B Z F Nl -2,N 4 -2(P)F* N2N3 ( P )+ 

A n a \ D%D%A%F NlNt (p)F* N2 _ hN3 _ 1 (p)+ 
D n a \ A n a \ A n a lDZF Nl -i,N 4 -i(p)F* N2 <N3 (p)+ 
BZCZCZBZF Nl - 2 ,N i -2(p)F* N2 _ 1>N3 _ 1 (p)+ 
CZBZBZCZF Nl - 1 , N4 - 1 (p)F* N2 _ 2tN3 _ 2 (p)+ 
BZDZDZBZF Nl - 2>Ni - 2 (p)F* N2 _ UN3 _ 1 (p)+ 
DZBZBZDZF^-^-^F^^p)]. 

Here N — 1 = (n — 1, a£cr), and -F/v'iv(q) is related to the 
Fourier transform of r} nq (r) in Eq. (J3j), 



n! / (-q y + zg^Rg 



(rc')I 



V2 



Appendix B: The exchange energy cost of the 
l z — ±1 excitations 

The transitions from n = to (n = l,a = 2) orbitals 
involve a large energy cost due to exchange with an in- 
finity of filled Landau levels. If an electron is promoted 
from an n = orbital to an n = 1 Landau orbital in either 
the ferromagnetic state |f.m.) or the valley-polarized one 
|v.p.) (l z = 1 excitation), the exchange energy with an 
infinite number of completely filled Landau bands (n > 2, 
a = 2,3, all combinations of £ = ±1 and a =t>i) con- 
tributes to the energy shift A®* a=2 ^(ofo-')? an analo- 
gous statement holds for the l z = —1 transition. This 
number may be formally infinite because we pick up a 
contribution from very low energies where the Hamilto- 
nian in Eq. (pQ) is no longer valid. This calls for some 
renormalization procedure. Exploiting the particle-hole 
symmetry of undoped bilayer graphene, however, this 
turns out to be a very simple one. Consider the fer- 
romagnetic state |f.m.) for concreteness. The exchange 
costs of the l z = 1 and l z = — 1 transitions are 



a ex _ /vex 



a ex _ _ a ex 

AA (0^cr)(l,a=2, ( e / cr / ) — ^oo 



dq 
dq 



C(p). 



(2t) : 



t-Iq; 



=2,la=2/ \ 
=2,la=2 IPJ' 



respectively. As the state obtained by moving an electron 
from the (0£ |) band to the (l,a = 2,£ t) band is the 
particle- hole conjugate of the state that is obtained by 
moving a hole from the (1, a = 2, £ t) band to the (0£ ^ 
) band, and the number of identical spin and identical 
valley bands is the same in the two cases, 



T la=2,la=2 
(27r )2 h a =2,l a =2 



(P). 



In the special case of the two band model this con- 
sideration yields = —^-f 1 A „f_ — ; here the infi- 



47reoe r ^B ' 



nite summation in A^ can actually be performed and 
yields an identical result. In the four-band model, on 
the other hand, the corresponding value ranges decreases 



from 



-0.16e^ 



to 



in the interval from B = 1 



to B = 40 T, and shows a tiny interlayer energy and 
valley £ dependence. 
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